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Abstract 

In the Randall-Sundrum (RS) brane-world model a singular delta-function source is matched by 
the second derivative of the warp factor. So one should take possible curvature corrections in the 
effective action of the RS models in a Gauss-Bonnet (GB) form. We present a linearized treatment 
of gravity in the RS brane-world with the Gauss-Bonnet modifications to Einstein gravity. We give 
explicit expressions for the Neumann propagator in arbitrary D dimensions and show that a bulk GB 
term gives, along with a tower of Kaluza-Klein modes in the bulk, a massless graviton on the brane, as 
in the standard RS model. Moreover, a non-trivial GB coupling can allow a new branch of solutions 
with finite Planck scale and no naked bulk singularity, which might be useful to avoid some of the 
previously known “no-go theorems” for RS brane-world compactifications. 
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1 Introduction 


In recent years, there has been considerable interest in the dynamics of brane interactions and brane- 
worlds with warped extra dimensions. The interest was motivated by the ideas coming out of the 
Randall-Sundrum (RS) brane-world scenario with a warped fifth dimension [1, 2]. The RS brane- 
world scenario promptly received many generalization in higher dimensions [3, 4, 5, 6 , 7] which have 
subsequently attracted much interest in gravity and cosmology research. A necessary ingredient of 
the brane-world models is that the space-time metric contains a warp factor that depends 

non-trivially on the extra dimension(s), and, as a result, a brane model of two 3-branes with opposite 
brane tensions [1], known as RSI model, may provide a geometrical resolution of hierarchy problem. 
This proposal was made more concrete in the scenario pioneered further by Randall and Sundrum, 
known as the RS single brane (or RS2) model [2], where the fifth dimension is non-compact. The latter 
model, viewed as an alternative to Kaluza-Klein compactification [2], gives an illustrative example of 
localized gravity on a singular 3-brane. 

To localize gravity to the RS brane, one considers a bulk space-time with negative cosmological 
constant. The low energy effective action in five space-time dimensions is taken to be 

S= cPxVW\(--2a)+J2 vW (£™ - ^. (1) 

Jb / i JdB 

where K 5 = 167 rG( 5 ) = ^^^1 with M(- 5 ) being the five-dimensional mass (energy) scale, and h is the 
determinant of the induced metric on the 3-brane. In a single brane setup, the brane tension T is 
positive (T >0), and the anti-de Sitter length scale I is set by a relation (^ = — 6 M^^g^/A. One of the 
interesting features of the RS action (1) is the presence of four-dimensional gravity as the zero-mode 
spectrum of a five-dimensional theory. One also finds correct momentum and tensor structures for the 
graviton propagator due to the “brane-bending” mechanism [4, 7, 8 , 9]. 

The action (1) describes the dynamics of a background metric field for sufficiently weak curvatures 
and sufficiently long distances. To further explore the general properties of brane-worlds, it is more 
natural to consider the leading order curvature corrections as predicted by string (bulk) theory. With 
warped space-time metrics in the bulk, however, one should take the curvature corrections in the 
RS action to be no larger than the second derivatives of the metric. Thus one can introduce the 
higher order corrections only in a special Gauss-Bonnet (GB) combination. The GB term arises as 
ol corrections in bosonic string theory [10, 11] and in heterotic M-theory scenario of Horova-Witten 
type [12]. Such corrections might be crucial in space-time dimensions Z) > 5, in particular, when the 
brane-world scenario is viewed as a low energy limit of string/M theory. 

There are now growing interest in the RS brane-world models modified by higher derivative correc¬ 
tions. Such corrections in a Gauss-Bonnet form, for constant dilaton fields, had been considered earlier 
in Refs. [13, 14], see also Ref. [15] which give some realizations of brane-world inflation due to quantum 
correction. It is learnt that brane-world configurations with a GB term and several co-dimensions one 
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branes can induce brane junctions of non-trivial topology [16, 17]. Furthermore, the presence of a 
GB term coupled to a bulk scalar in the effective action leads to interesting physics in a variety of 
context, ranging from gravity localization [18, 19, 20, 21, 22, 23, 24, 25, 26] to FRW type cosmology 
on a brane [27, 28, 29]. It is worth noticing that a five-dimensional brane-world model with a GB 
term reproduces all essential properties of the original RS models, including the tensor structure of a 
massless graviton [9]. In this paper, we extend the work in Ref. [9] to the D dimensional space-times, 
and also examine the general properties of brane-world solutions without and with a bulk scalar. 

A fine tuned relation between the bulk cosmological constant and the brane tension is required 
in the RS model to maintain flatness of the 3-brane [30]. This problem has been known as ” no-go 
theorem” [30, 31, 32] for non-singular RS or de Sitter compactifications based in Einstein gravity. 
It is therefore of interest to know whether these arguments can be changed to cover alternatives in 
the brane-world actions. With a GB correction to Einstein gravity, fine tuning is required in the RS 
models to get a singularity free solution with the finite Planck scale [33]. Nevertheless, in the presence 
of such interaction, there exists a new branch of the solutions, for which it might be possible that 
the only bulk singularity occurs when the warp factor vanishes at the anti-de Sitter horizon, 

z —> oo. Presumably, a GB term might smooth out the bulk singularities and hence avoid some of the 
previously given no-go arguments for the RS compactifications. 

The rest of the paper is organized as follows. In Sec. 2, we introduce the brane-world action and 
find intersecting brane backgrounds with more than one uncompactified extra dimensions. We give 
the basic expressions of the linearized Einstein equations with a warped metric, and generalize results 
of the RS2 model when there are two transverse directions. We then discuss in Sec. 3 some interesting 
features of the RS solutions modified by a GB term. In Sec. 4, we find the Neumann propagators in 
D dimensions, where the background is described as a 3-brane embedded in AdS space of dimensions 
D > 5. In Sec. 5 we analyze the energy conditions and compare the RS solutions in five dimensions 
with those arising due to a non-trivial GB coupling. In Sec. 6 we give some insights on the nature 
of singularities or no-go theorem for a class of brane-world gravity models coupled to a bulk scalar 
field and a Gauss-Bonnet self-interaction term. Section 7 contains discussion and outlooks. In the 
Appendices, we give some useful derivations of the linearized equations for a class of higher derivative 
gravity in brane backgrounds. 

2 Gravity in Brane Backgrounds 

We shall begin with the following H-dimensional gravitational action 
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where dB represents the {D — l)-dimensional boundary. The gravitational coupling kd = = 

with M and Gd being, respectively, the D dimensional mass scale and Newton constant. The 
indices {p, q, ■ ■ ■, s) = (0,1,2, 3, • • • , D — 1), A is a Z) dimensional bulk cosmological term, {k = 
1,2, ■ ■ ■, (D — d)) represent the vacuum energy (brane tension) of the (D — 2)-branes, and T is the d 
dimensional brane tension at a common brane junction. For practical purposes we shall take d = 4, so 
a takes a mass dimension M^~^. The second action in (2) is the effective action for {D — 2) branes, 
while the last term introduced at the common intersection of higher dimensional branes characterizes 
a four-dimensional brane action. This term will be in effect only if D > 6, because in D = 5 the 
second action in (2) has already represented the sum of 3-brane actions. For vacuum branes, one has 
= 0, and also — 0) since the matter degrees are supposed to be confined on the branes. 


2.1 Choice of background and RS tunings 

Let us consider a smooth version of the multidimensional patched AdS space, introduced in Ref. [3] 
to study intersecting brane-world models, with metric 

dx^ dx'^ + gij dz^'dz^) , (3) 


where p, iz = 0,1, ■ ■ ■ {d — 1) with N = {D — d) being the number of extra spatial dimensions. The 
Einstein field equations {a = 0) give a solution A{zi) = log(X)£=7^ kil/f’+l)) with i being the curvature 
radius of AdS space, by satisfying the RS relations 

{i)-l)(D-2)(D-4)Af»-2 , 4(D-2)A4°-2 

A =-^2-. f'^■ = - i -■ W 

One has T = 0 for a = 0. The four-dimensional brane tension T at the common intersection of two 
4-branes is non-zero for a > 0 [16, 17]. The two expressions in (4) imply the RS fine-tuned relation 


A 


2 

k 


32(11-2) 

{D-1){D-A) 


A. 


(5) 


Therefore, since A| > 0, the bulk space-time is anti-de Sitter (A < 0). For D = 5, can be replaced 
by T, and hence A^ defines a 3-brane tension in an AdS^ space. 


2.2 Linearized Einstein gravity 

Next we consider a linearized theory. With warped space-time metric (3) in the bulk, the linearized 
R-dimensional Einstein field equations take the following form: 

_la2^ -L-n QCQ D — 2 \\f)U P) h — r) h 

2 ‘^p(J ' 2 2 Gf Gfpitqra t 

[2d,nd^A -{D- 3)d^A d^Aj hp^ 

[2d-d”A -{D- 3)d^A d^A] rip, 

= -kd Ahppe-^^^^^, ( 6 ) 
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where we have used a harmonic gauge d^^hpq = ^ /if. Since the warp factor A{z) is only a function of 
the extra spatial coordinates, the indices m and n take face values from zi to ZD-d- When d = 4, the 
number of extra dimensions N is {D — 4). The components of ( 6 ) largely simplify in using the 
RS fine-tuned relation (4). The second term in ( 6 ) would not show up in using transverse-traceless 
gauge V^hpq = 0 = hp, instead of the harmonic gauge. 

It is more convenient to perturb the metric background (3) in the following form 

ds^ = {{VtJ.i' + hpu) dx^ dx^ -\- dz^dz^) . (7) 


However, in this gauge, other than the standard graviton hpu, there may be some additional grav¬ 
itational degrees of freedom coming from hpm and hmm- With more than one (conformal) extra 
dimensions there is a general problem for diagonalizing the linearized fluctuations of all metric fields, 
which is computationally difficult. To this end, since a gravitational potential on the brane is mediated 
by effective 4i4-gravitons, we find reasonable to analyze only the linearized field equations for graviton 
with the gauge /i^ ^ = 0, d^h\p = 0. One may use an approximation [3, 6 ], where the tensor modes 
either decouple from the off-diagonal components of the perturbed equations like hpm and diagonal 
components hmrm or only the nonzero components of the fluctuations are hp^. We comment upon the 
decoupling of the scalar modes of the metric perturbations in the discussion section. 

In our analysis we follow the background subtraction technique introduced in Ref. [ 6 ], where one 
considers vacuum branes and subtracts out the background fields from the variations of Gab and Hab 
(the second order Lovelock tensor). For a gravity action of the form (2), one may define dTab = Tachl, 
where Tab = Gab is taken about the background. This approach extends in an obvious way to 
developing the perturbative expansion to higher order curvature terms. An additional benefit of this 
approach is that the RS fine tunings of the previous subsection will be only implicit. Furthermore, 
it is reasonable to impose the gauge dph^^ = 0, hp = 0. With these approximations, the Einstein 
equations 5Gab = i^ndTab linear in hpy take a remarkably simple form [ 6 ] 


- dad^hp, + {D-2) d^A dmhp, = 0. 


( 8 ) 


We define hpy = and ^4 hpy = m^hp^, where hpy = epue^P'^ip{z), Cpu is the polariza¬ 

tion tensor, and arrive at the following analog non-relativistic Schrodinger equation 



D{D - 2) 


D-i 

^ sgn{zif 
i=l 


P-2) 

i\zi - 1 - i\) 





m^'4){zi). 


(9) 


We are considering here the case where warp factor ig a function of all transverse coordinates Zi, 

where i = 1, 2, • • •, (74 — 4), so that Xlf sgn{zi)^ = (74 — 4) other than at the origin (brane-junction) 
in the transverse space. In order to normalize four-dimensional metric at the origin, one has to assume 
that A(0) = 0. For 74 = 5, Eq. (9) gives the RS one-dimensional Schrodinger equation, and this was 
extensively studied, for example, in Refs. [4, 7, 6 ]. So we shall be interested only in the 74 = 6 case. 
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Define a set of new coordinates x± = l^il ± \z2\, so that the bulk part of (9) takes the form 


- dz - dz + 


6 


"+ {\x+\ + \x.\+^)\ 


i>{x-,x+) = l-m^^{x-,x+), 


( 10 ) 


where ip{x-,x+) = V’(- 2 i,-Z 2 )- We can separate this equation into 


-dl_(p{x-) = m‘i(p{x-), and (p{x+) = ml (p{x+), ( 11 ) 


where we have defined = 2(m^ +m^), and V’m = ^m- (x-) x ^m+ {x+)- Since > 0, m?_ and ml 
each can have only positive eigen energy. The solution for each continuum wavefunction is a linear 
combination of Bessel functions 


(pm_{x-) = Ni{m-)\/\x-\ (sin(m_x_) + Amcos(m_x_)) . ( 12 ) 

iPm+{x+) = N2{m+) ^/{\x+\ +€) [J 5/2 (m+(|x+|+£))+ 5 mT 5/2 (m+(|x+|+£))] . (13) 

These solutions must satisfy the Neumann type boundary conditions at the brane junction x± = 0: 

[x±-^±^^ + ^i^{x±) = ^- (14) 


The coefficients Am and Bm can easily be read off 

Am = cot(m_x_), Bm = 


43/2 (?^+^) 

•^3/2 (" 1 +^) ’ 


(15) 


Thus V'm+(0) ~ (m+£)“ where a + 1/2 = 5/2. There are 4 (= 2^ ^) sectors in the mass spectrum 
spanned by (i) ipo{x-) (po{x+), (ii) ipm_{x-) ipQ{x+), (hi) ipQ{x-)ifm+{x+), and (iv) ipm_{x-)ifm+{x+). 
The state (i) gives the four-dimensional graviton, and the set of continuum states (if) or {in), which 
is localized in x_ or x+ direction, contributes as an integral over the single eigenvalue m_ or [ 6 ]. 
The set of continuum states (iv) contributes to Newton’s law. For two point masses mi and m 2 placed 
at a distance |x — x'| = r on the four-dimensional brane intersection, the Newtonian potential is 


—A4l/(r) ~ 


mi m 2 1 




V^m_(0)|^ 


V^m+( 0 )| , 


(16) 


where m = y m_ -|- ml . A qualitative feature of this potential can be known by specializing to the 
case where m_ would extend down to mo = 0, and V’m_(0) = 1. Hence 

_A.y(,.) = 0W=i^..(l)/ (17) 


where l/G^^ ~ Cl is a number of order one. For large D and i < r, the Kaluza-Klein mode 

corrections {(./r)^~^ are more suppressed, so gravity becomes weaker as the number of transverse 
directions increases. Similar results are known for gravity localized on a four-dimensional string-like 
defect in D = 6 [34], a difference here is that now we have two (non-compact) transverse directions, 
and this is itself more interesting in the scenario of Refs. [35] and [3]. 
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3 Corrections to Einstein Gravity 


Next we investigate the brane-world solutions which occur due to a non-trivial Gauss-Bonnet coupling 
a. The exact metric solution for the modihed Einstein equations is given by A(z) = log(^^ \ zk\/L + 1), 
where L is the AdS curvature radius which has a contribution from the GB coupling. We should note 
that L? is defined by the bulk solution [17] 

1 1 

l2 “ 2{D-4.Y{D-?,)aKD 


l± \ 1 + 


8{D — 3) {D — 4) aA k 


D 


{D-l){D-2) 


(18) 


The space-time metric is therefore 


dsf) — 


1 


{rj^ydx^dx'" -f ^ {dzjf^ . 
1=1 


(19) 


One may view the above metric background as a four-dimensional intersection of higher dimensional 
branes or as a 3-brane embedded in space of dimensions D > 5. One has = 2{D — 4)^(Z1 — 3) a kb 
when the two branches of the bulk solution (18) coincide. The RS type background relations 


(I 1 - 1 )(I 1 - 2 )(Z 1 - 4 )( 1 - 7 / 2 ) ^ _ 4(74 - 2) (1 - 7 / 3 ) 

2L‘^kd ’ ^ Lkd 


( 20 ) 


may be assumed implicitly for 7 < 1. Here 7 = 2(74 — 3)(74 — and a is the usual 

Gauss-Bonnet coupling, so 7 = 0 when a = 0. In this limit, one finds from (20) the RS fine-tuned 
relations. We may define a = a\, with ai being a dimensionless GB coupling. If we define 

I = L/^/{D — 4)(1 — 7 / 2 ), then A takes the usual form A = — {D — 1)(74 — 2)/(2kd l'^)- 

Let us include a matter source on the brane, and study linearized equations for the effective 
four-dimensional gravitational fluctuations. We consider only the scalar wave equation for each of the 
components in the background (19). The linearized field equations for read 


1 - 


2(74 — 3)(74 — 4)^ ano 


7.2 


,2 q 2 I 2 ) 


-di-dt+ ^ 


sgn{zi)dzAh fii/z) 


2{D — A) aK£) xt 

H-r-/ 

i¥=j 


L 


-dl- sgn{zk)d,^ h^^ix, z) = -kd ■ ( 21 ) 


One of the important features of the modified solutions due a non-trivial a is that a Gauss-Bonnet 
term contributes with a delta-function term to the linearized equations, which implies a non-trivial 
topology at the brane junction, like that a non-zero four-dimensional brane tension at the common 
brane-intersection with several co-dimension one branes. In Eq. (21), terms involving Dirac delta 
functions vanish for \z\ > 0. Therefore, for z ^ 0, one has a RS type bulk equation but multiplied by 
the factor (1 —2(74 —3)(74 —4)^ okb/L'^) = 1~7- This apparently implies that the effect of a GB term 
may be removed by redefining the Einstein constant. But this is not the case, rather a non-trivial GB 
coupling a alters the physics of the brane-world by modifying the Neumann propagator in the bulk. 
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If one allows the parameters to take values such that the first square bracket term in (21) vanishes, 
the solutions become solitonic. This is because what survives after this setting precisely gives a RS 
equation multiplying with an additional delta function term, but in one co-dimension lower. As in 
the RS solution, which is given by 7 = 0, there always exist brane-world solutions for 0 < 7 < 1. So 
the physical relevance of the GB coupling a is two fold. If one has solutions satisfying 0 < 7 < 1, 
such solutions will be fairly similar to the RS solutions except with some corrections, like in the 
Newtonian potential. But if one is allowed to take 7 ~ 1, one finds a new branch of solution with 
finite effective gravitational constant without finite distance bulk singularity. In the next section, we 
shall be interested in the 7 < 1 solution. 


4 Green Function in D Dimensions 


In this section we analyze the Green functions by expressing a L)-dimensional propagator Q£)[x, z] x', z' 
in terms of the Fourier modes such that 


Gd{x,z-x',z') = J 

where the Fourier components gp{z,z') (= h{ z, z') ~ 'G{z,z')) satisfy 

(1 - 7) d,)Gp{z, z') = e(^-2)^G) - z'). 


( 22 ) 


(23) 


When one of the arguments is at z' = L, the Neumann propagator is calculated in Appendix D to be 


Gd{x,z-,x' ,L) = (1 - 7 ) 


-1 


l) J (^ 


\2u 


Jp(x-x') 


O 1 


Hi^\qz) 


nHi^},iqL)+xiqL)Hi^\qL)i 


(24) 


where u = {D — l)/2, x = — 3)(1 — 7 )), and = J + iY is the first Hankel function and 

q"^ = —p^Pp = m?. For both the arguments of the propagator at z, z' = L (on the brane), using the 
Bessel recursion relations, one can find a two point correlator of the effective theory. For example, 
when D = 5, a two-point correlator is 

2 I (1-7)<^((?T) 


{(j){p)(j){-p)) = 


(25) 


(1 + 7 ) [q^L g (1 + 7) //W (qL) - 7 (qL) 

For arbitrary D, by using Bessel expansions in Eq. (24), we obtain the scalar Neumann propagator 

■{D-3) 


Gd{x,L;x' ,L) = (1 -^ 7 ) 


-1 


L 


■ Gn-iix, x') + Gkk{x, x') 


where 


Gd- 


Xx,x') = j 


d^~^p 

(27r)^-i 


(26) 


(27) 
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Gkk{x,x') ~ -j 


jD—1 


(27r; 


P ^ip(x—x') 


D-l 


{D-^ + {D- 6 ) 7 ) 1 


2{D-b) (1 + 7) 


qL 


1-7 {qL)^-^ 

1+7 qCi V 2 


( 28 ) 


where Ci is a dimension dependent number (see App. D). For 7 = 0 , the above results reproduce 
the scalar Neumann propagator found in Ref. [7], after a substitution d = D — 1, and the results in 
Ref. [9] with D = 5. The first term in (26) is the standard propagator of a massless scalar field. 

The long distance behavior, r >> L, of the propagator is governed by a small q behavior of the 
Fourier modes. Thus, for qL « 1, a leading order contribution to the propagator comes from the 
logarithm. For D = 5, and \x — x'\ » L, qL « 1, we find, to the leading order in g. 


Qa{x,x') = j 


d^p 


Jp{x-x') 


1 


OC 


\X — X‘ 


/|2 


(29) 


Thus 1 / 4 (x, x') is just the ordinary massless scalar propagator in four-dimensions. Eq. (26) implies that 
a four-dimensional massless graviton propagator is contained in the full five-dimensional propagator. 
This is one of the plausible results of the RS2 model supplemented by a GB term. By the same token, 
for the KK modes, one has 

l-7rD-4 /■ 


Qkk{x,x') ~ 


OC 


1 + 7 
1-7 


L 


/ 


(27r 


p 


^ip(x-x’) ^(D-5) in(gL/2) 


LD-a 


(30) 


1 + 7 |x — 

For 1 > 7 > 0, we expect some correction in Newton’s law. At large distance scales on the brane, 
|x — x^l = r » L, the contribution of the KK modes (30) is very small compared to the zero mode 
contribution (26). Like the 7 = 0 solution [7], for H > 5 and odd D, a leading behavior of Gkk goes 
like ~ ^ 6 - 2 D^ When D is even, there are no logarithmic terms and a leading behavior of Gkk{x, x') goes 
like ~ Lr^~^, which is further more suppressed compared to the zero-mode contribution. Therefore, 
the effect of gravity becomes weaker when the number of extra (transverse) dimensions grows. One 
also notes that the Newtonian potential for a point source of mass m, and H = 5, is 


K(r) ~ — 


G^m 


1 + 1 - 


27 

1 + 7 


1 ^ 
2 


(31) 


In this formula, the contribution of the brane-bending mode is not included, which might change the 
factor 1/2 into 2/3 [4, 9]. It is seen that the GB term contributes to Newtonian potential with the 
opposite sign to that of the Einstein term or scalar curvature. One can expect a small correction to 
Newton’ law with 7 < 1 and r > L, but such a correction is almost trivial for r » L or/and 7 ~ 1. 


5 Nonperturbative Solutions with a GB Term 


A non-singular Minkowski or de Sitter brane-world compactification could be difficult for Einstein’s 
theory [31, 32]. So one is lead naturally to hope that adding higher derivative corrections to the 



brane-world action might improve this situation. As a convenient choice, we shall begin with the 
Lagrangian of gravity, including a GB term and a bulk scalar field, in the form 

S = Jd^x ^ {d^f - 2V{^) + •••). (32) 

The dots in (32) represent higher order terms in the scalar field ((/?) and some other fields which 
are not turned on. We assume that V((f) is non-positive (bulk) cosmological potential, and (p can 
only depend on extra coordinates, p = p{z), as dictated by the Poincare invariance on a brane. 
Though we will not restrict the action (32) to a particular string theory background, one may think 
about it like including a' corrections in type I string theories. In propagator correction-free Gauss- 
Bonnet scheme [11], the /^^-corrections can have a dependence on a scalar field, a ^ ^ Aoe“™'‘^, with 
m = 1/\/D — 2, Ao = 1/4 (1/8) for bosonic (heterotic) string. For a constant bulk scalar, since V{(p) 
takes a bare value Aq, this parameterization is not so important for our analysis. 

We can write the metric background in the form 

= n‘^{z) (dx^ -|- gmndz"^dz"') . (33) 

Here 11(2;) = dx^dx'^ with g being the d-dimensional metric, which can be Minkowski, 

de-Sitter or anti-de Sitter space. We study only the case of H-dimensional warped geometry compact- 
ified to d-dimensional RS type space-time {g^^ = rjfjiu)- 

The equations of motion in H-dimensions take the form 

Rabid) — ^ 9abTc^ T Rab — bC£) T^b , (3^) 

where the contribution to the stress energy of a massless scalar field and a bulk cosmological potential 
reads 

Tab = dapdbP - -gabidpf -V{p)gab , (35) 

and, the contribution of the GB interaction term reads 

Rab = _2 dabRcB “ ^Q; (^RRab “ ‘^RacbdR^'^ + RacdeRb “ 2Ra '^Rbc^ ■ (36) 

5.1 Energy conditions with a Gauss-Bonnet term 

Let us take the {gf ) components of the D dimensional Ricci tensor 

= d?,..(d)-d/.. (V^ log L! +(11-2) (V log 11)2) _ (37) 

In the RS compactification, d-dimensional space is a Minkowski space, so R{g) = 0. In our metric 
convention (mostly positive), for a = 0, the ”c-theorem” of Ref. [31] reads 

-(d-l)e2^W(A" + A'2) = R\-Rl = n,+i{Tl-Tl) 

= -Kd+i < 0 , (38) 
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where V(99)’s contribution to r/ and cancels out. This is the obvious positive energy condition 
Tab ^ 0, where is any arbitrary future-directed timelike or null vector. Since 

= d{d-2){d-3)a(2A''+ A'‘^^ (39) 

there arises an identical condition from 

Kd+i (T/ - r/) = 2e{d - (^A" + A'^^ > 0, (40) 

where e = {d — 2){d — 2>)aKd+i- The condition A” + A'"^ > 0 may be equivalent to the weak energy 
condition (WEC) —p + \pi\ > 0, with p and p being the energy density and pressure. For a > 0, the 
strong energy condition (^Tab — dabT^ > 0 is promoted to Tab > 0. First consider 

(^Tab - 9ab = (^daP> d^ip - V(t) ffab^ 

= {d,pf eiz - v{ip) (e%+ eiz) > 0. (41) 

So the strong energy condition is not violated if is a null vector and V(p) < 0. For a future-directed 
null vector, = 0; Eq. (41) implies that — {dzp)^ > 0; so can be 

timelike (^^^^ = —1) or null-like (^^^^ = 0). This is consistent with the condition Rab^^'C^ ^ 0^ 

-(d - 1) + A"'^ > 0 . (42) 

We would like to see whether the second piece in (34), i. e., the contribution to stress energy from 
the GB term, violates the strong energy condition. For a null vector the condition Tab^°'^^ > 0, 
using (39), yields 

(-2a {d - l)(d - 2){d - 3) (yl'^ + yl") > 0 • (43) 

Thus for all non-spacelike i. e., < 0, the strong energy condition (SFC) is intact provided 

that the WFC yl'^ -|- yl'^ > 0 holds. This may not be the case for any other combination of Riemann 
tensors or terms Here, validity of the SFC basically says that brane gravity is attractive. 

5.2 Randall-Sundrum limit 

Let us introduce a brane action for the RS singular 3-brane with a positive brane tension cj > 0: 

Sbrane = 2 J d*X ^/\^\ (-fj) . (44) 

^We acknowledge fruitful correspondences with C. Nunez about the no-go theorem and energy conditions with a GB 
term that prompted us to add the above explanation. 
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For a constant scalar, the field equations following from (34) in D = 5, including a contribution from 
(44), simplify to 


3 (^A'^ + A") 

12 ^ 4 '^ 


3 (^A'^ + (2s + a K 5 Kz) 

- K 5 Ao + 6 (2e A'^ , 


(45) 

(46) 


where e = 2aK^. For e = 0, the bulk equation A" + = 0 may imply that (= C) 

is a constant. However, C is already fixed by Eq. (46) such that C = — where £ is a 

constant with dimension of length. Solving the equation A" + A'^ = 0 amounts to selecting a solution 
A'{z) = (l^l + z*)~^ of the unperturbed Einstein theory, but z* is undetermined by the bulk equations. 
One may fix z* using the continuity condition and normalization condition, such that A'(0) = l/£, 
and arrive at [2] 


a = 




Ao = - 


12 


K5 


£2 


(47) 


As is well known one has some undesired features with A" + A'^ = 0 as an initial condition. For 
example, there is a finite-distance bulk singularity for z* < 0 at z = z*. So we will now investigate the 
case of interest, e > 0. In this case, the bulk equation can be satisfied even by selecting 


2ee2^(^)A'^ = 1. 


(48) 


It is obvious that this freedom is not there with e = 0. Since e > 0, the metric solution is 


0 ^( 2 ) = 


1 




dz 


„-A(2) ^ 




+ v/^’ 


(49) 


where we have normalized the solution A(0) = 0, such that e takes a value 1 at 2 : = 0. We find 
it convenient to define the length scale = /, so A^' -|- A'^ = The condition (48) therefore. 


unsurprisingly, also the solves the bulk equation A" + A'^ = 0. It is obvious that A" -|- A^^ = j on 
the brane (2 = 0 ), as it should be in order to keep the WEC A" -|- A'^ > 0 intact. The solutions with 
e > 0 have all essential properties of the RS solutions, such as converges as 2 : ^ ± 00 , and there 

is no any finite distance bulk singularity. However, there is now a new length scale in the problem, I, 
and a common feature of these new solutions is that they are not analytic in the coefficient a of the 
Gauss-Bonnet interaction, so e > 0 is a physical requirement. Remarkably, the RS solution with e = 0 
corresponds to a limit of these solutions where the singularities are pushed to infinity, 2 ^ ± 00 , so 
are harmless as they might have interpretations in field theory as (ultraviolet) energy cut-off scale. 

The boundary condition relates I to the brane tension a. For a brane at 2 = 0 with positive 
tension, a > 0, Eq. (45) implies that 


65{z) 
I + \z\ 


I - ^ sgn{zf ] = aK 5 5{z). 


This after regularizing the (5-function, 5{z) sgn{z)‘^ = 5{z)/3, determines the brane tension 


a = 


1 4 

^5 I 


(50) 


(51) 


II 



The tension a is generally not fine-tuned because I is arbitrary. However, the scale I is used also to 
fix the bulk cosmological constant, via Eq. (46), 

= (“> 

There is a fine-tuning between a and Aq, which is required to maintain flatness of the 3-brane. This 
may be relaxed for more general solutions, like (anti-) de Sitter branes, with e > 0. 

A remark is in order. The bulk scale I is somehow fixed by the coefficient e, but this is not an 
unnatural choice from the view point of low energy effective string action, rather a common result 
in higher-curvature stringy gravity, see, for example. Refs. [14, 19, 20]. In the Einstein frame, the 
tree-level (bosonic) string action, with appropriate conformal weights for a dilaton (f), reads [11, 39] 

Sbuik = j d^x yj-go + Aoa^ e~ ^ _ 2 

- m , (53) 

where m = 4/(Zl — 2). Therefore, with (j) = (po = const and D < 10, the bulk cosmological term 
Aq oc —Ija' ~ — 1/e. The above action diverges at a' = 0, so one might require a' > 0, so e > 0, 
in order to get solutions which are free of singularities. In the context of lowest-order brane-world 
Einstein gravity, so a' = 0, there arise naked bulk singularities due to the fact that both dilaton and 
graviton field exhibit logarithmic singularities [40]. This problem can be resolved by considering the 
leading order a' corrections as recently shown in Ref. [39]. 

6 Brane-World No-Go Theorem 

Basically, with warped space-time metrics in the bulk, there are three different arguments given for 
the brane-world “no-go theorem” in the Einstein theory. They are 

• No singularity free solution with a finite Planck mass is possible without a fine tuning [30]. In other 
words, singularities in the self-tuned solutions are generic if gravity is to be localized. 

• It is impossible to have A'{z) approach a positive constant as z ^ 00 and a negative constant 

as z ^ — 00 . This monotonicity of is often called the brane-world c-theorem [31] 

• There are no non-singular Randall-Sundrum or de-Sitter compactifications where the only possible 

singularities occur when the warp factor goes to zero at the singularity [32]. 

The first argument above is mainly related to the fine tuning of the cosmological constant, which 
may complement the Hawking-Penrose singularity theorem [36], rather than being intrinsic to the 
Randall-Sundrum models It is interesting to know whether any of these arguments can be avoided. It 

^In terms of the y coordinate such that dy = one has A'{y) = A'{z), A"{y) = (^A"[z) + A'^(z)^. 

®In the RS brane-world context, this implies that the generic initial conditions, such as A"{z) + A'{z)^ — 0 for z > 0, 
lead to singular solutions of Einstein equations. 
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should be possible to circumvent some of these arguments by considering higher derivative corrections 
to the gravity equations [32], 

6.1 No-go theorem and possible avoidance 

Let us take the (fiu) components of (34) and contract the ^j^-indices to arrive at 
d {d-l)A'‘^-A" = 

-d{d-2)aKd+i 2A"-{d-2)A'^ (54) 

where 

+ = (55) 

Here we have used the fact that ip = ip{z), and assumed V{ip) < 0. 

First we briefly review the no-go theorem advocated in Ref. [32], With 0 = 0 and d = 4, one has 

2,A'‘^ - A" >0. (56) 

Integrating (56) over the compact internal space by parts one finds f dz\/g < 0. This 

may be satisfied only if one allows the equality sign in (56) and is a constant. The condition 

Roo = Too — 75 ^ ffoo Tc = 3H'^ — = 0, or Rqo = T = 0, may not generate gravitational fields in the 

bulk space-time [37]. Moreover, in the RS compactification should not be a constant, rather this 

factor is essential to explain the warped nature of a bulk geometry and hence a RS compactification. 

We will now investigate the case of interest, a > 0, so e > 0. For d = 4, the four-dimensional parts 
of the curvatures read 

R{x) = (x) -h 4 (2A" - {N + 2)H'^) 

‘^{N + l)A"{z) 

- 2N{N + l)A'^{z)\ + • • • , (57) 

where N = D — A. In the following we restrict our attention to D = 5. The four-dimensional Planck 
scale is therefore [17] 

/ OO P . . -| 

dz 1 + 40^5 (2y4" - 

/ OO , p -| _j_ QQ 

(l +40^5 +MJ) 80^5 . (58) 

R is obvious that Mp^ is finite for A{z) = In (1 -|- | 2 ;|//). From Eq. (54) one has 

4 (SA'"^ - H") = V{ip) - 8e (^A" - . (59) 
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Thus, with a > 0, the inequality (56) may hold even in reverse order without violating the positive 
energy condition, and the warp factor e~ need not be a constant. In particular, when the weak 
energy condition A^' > —A!'^ saturates, we find 

16 (^1 - e ^ . (60) 

This can be satisfied when I^((/?) < 0, > 0, and < 1. We have shown that in the presence 

of a GB term there are some improvements over the no-go theorem. 

6.2 Solutions with a bulk scalar 

For a non-constant scalar field, the bulk equations of motion following from (34) simplify to 


/2 

l^d+l 

= 2(d - 1) (^A'^ + A"^ (^1 - 

(61) 

2Kd+i V ((p) 

= -e^^^^\d-l)\{d-l)A'^ - A"' 



+ {d- [(d - 2) A'^ - 2A"'\ , 

(62) 


where e = {d — 2){d — 3)aKd+i- We may express these equations in the y coordinate such that 
dy = e~^^^'>dz. For d = 4, and (p = ip{y), = A'{y) = W{p), dW{ip)/dip = W^p, Eqs. (61), 

(62) take the following form; 

( 3 \ 3 1 

-^-Ao (l-2eW^(^)f + Ao, Aq =^ , (63) 

(p'(y) = (1 - 2£W^((p)) . (64) 

K5 

The above (super) potential (63), named due to its resemblance with supergravity solution, for the 
e = 0 case, was analyzed in Ref. [38]. In the following we shall be interested in analyzing bulk solutions 
with e > 0, without retaining explicit form of A(y) (see Refs. [23, 25] for some relevant discussions). 

Some non-singular solutions are found in Ref. [39] by satisfying — 22.2 < aAg < —5/12, in the 
units K5 = 1. So in the bulk it might be possible to take V((p) = Aq and (l — 2elT^ ((/?)) = 0. 
In this limit, a domain wall solution smoothly interpolates between two anti-de Sitter minima (p± 
of the potential V((p) as y ^ ± 00 . Because p'{y) vanishes in the bulk, V((p) takes a bare value 
^(v^±) = ^ 0 - The scalar field and warp factor simply become ip = ipo and A{y) = ± y/— Aq \y\. One 
notes that 1 — 2sWp / 0 on the brane, rather this becomes 1 — 2elT^/3 > 0 due to an essential 
(5-function regularization. There is no any bulk singularity for the (super) potential of the linear form 
W{p) = Ai (/? -I- A 2 [14]. 

It might be desirable to know what would happen if 2elT^(e) 7 ^ 1 in the bulk? Then, since 

^'{y) = ^'(.y) (1 “ 2eW‘^{p)) , (65) 
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the A"{y) > 0 condition holds only if 2£W‘^{^p) < 1. In this case, there is no improvement over the 

c-theorem. As we expect that A"{y / 0) = 0, W^p should vanish in the bulk, and hence W{ip) becomes 

a constant function of scalar field, which implies that a bulk scalar is not dynamical. 

In the e = 0 case, the condition 

A"{y) = -^ > 0 ( 66 ) 

Kb 

is used to prove a c-theorem in Ref. [38]. With e = 0, the metric solution reads A{y) ~ \y\/^, where 
i = y^ 6 /(—A). In general, a domain-wall solution should interpolate between y —oo (infrared 
region) where A{y) is linear, and y ^ oo (ultraviolet region or AdS horizon) where A{y) is again 
linear. But the condition A"{y) > 0 rules out here the second possibility. So it is not possible to 
have A'{y) approach a positive constant as y ^ -|-oo and a negative constant as y — oo. This 
monotonicity of A'{y) has been known as brane-world “c-theorem” [31, 38]. However, with e > 0, and 
2sA'‘^{y) = I as a bulk solution (c.f. [46]), warp factor gets the both signs 

A{y) = . (67) 

Then it might be possible that A'{y) approaches a positive (negative) constant as y ^ -|-oo (— oo). 
Therefore, the brane-world “c-theorem” of [31] may not be available to the a > 0 case. 

7 Discussion and Outlook 

If the Randall-Sundrum models are to be the low energy limits of some fundamental (string) theory or 
yet-unknown theory of quantum gravity, it is likely that the gravity equations include higher curvature 
corrections such as a Gauss-Bonnet invariant. For a background of branes coupled to matter sources 
and a GB self-interaction term, we have presented some useful expressions for the Neumann propagator 
in arbitrary D dimensions and analyzed the structure of graviton interaction. It is shown that the 
RS model with a GB term in the bulk gives a massless graviton on the brane as in the standard RS 
model. Perhaps this is one of the most striking results of this paper. We have shown that for a small 
GB coupling a, so 7 < 1, the brane-world solutions are qualitatively similar to the RS solutions. We 
also pointed out a possibility that the Newton’s law is exact, other than that such a behavior one 
would expect at large distance along the brane, if one is allowed to take 7 in the order of unity. 

We have examined the general properties of the RS solutions coupled to a bulk scalar and a GB 
term and have found that fine-tuning is a generic feature of RS models. We analyzed the energy 
conditions with warped space-time metrics in the bulk and found that energy conditions are not 
violated by a GB term. Meanwhile, we observed a new branch of solution with finite Planck scale 
and no naked bulk singularity. More precisely, a bulk singularity for the e = 0 solution, which may 
complement the Hawking-Penrose singularity theorem [36], is pushed for e > 0 to the singularity at 
the anti-de Sitter horizon 2 ; ^ ± 00 , so is harmless as it might have held theory interpretation. The 
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new solutions with e > 0 are also found useful to avoid some of the previously known no-go arguments 
for the RS brane-world compactifications. 

In this paper, we have performed the calculations with an assumption that the scalar modes of the 
metric fluctuations decouple from the tensor modes. This is a possible one, because, at least for the 
flat RS branes, there are no delta function sources for the scalar and vector modes [17]. In the presence 
of several co-dimension one branes, such a decoupling certainly requires a moduli stabilization and it 
is possible that the physics responsible for this stabilization would modify the analysis for the tensor 
structure of the graviton propagator, which is not analyzed here due to this subtlety. This treatment 
might require more general assumptions that there is a bulk scalar coupled to brane gravity, and the 
branes are (anti-)de Sitter. We hope to return to this point in future publication. 
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Appendix A: Metric variations with quadratic curvature terms 


The starting brane-world action is the Lagrangian of gravity which has a general form in D dimensions, 
including quadratic-order curvatures. 


S = [ (fx^J-go 

— -2K+(aR^ + PRabR^’^ + -fRabcdR^^"^) 

Jm 

kd V / 

/•i’th brane 



JdM 


+ J dA^x^J-gd 

i-T). 


where a,b, - ■ ■ denote the L)-dimensional space-time indices. For D = A, the action (A.l) will be free 
of massive spin-2 ghost only if /I -|- dy = 0 [41]. For D > A, however, one requires that a = —fSjA = 7 , 
which is the Gauss-Bonnet relation. The graviton equations derived by varying the above action with 
respect to may be expressed in the following form 


\/—9D Gab + Hab + ^gab) 



D-4 




2 = 1 


T 

'2 



{zi,Z2,--,Zn=0} 




(A.2) 
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where Hab 


Hab = -^gabiaR^+/3RcdR^'' + lRcdefR^^^^) 

+2[aRRab + f^RachdR'^^ + l{RacdeRb ~ “^Ra ^Rbc + ‘^RacbdR'^'^)] 

— (2q; +/3 + 27)(VaVfei? — + (/3 + 47)V^ ^iiab — — • (^-3) 

The curvature derivatives vanish for a Gauss-Bonnet invariant (i. e., with da = —(] = 47 ). Then 
the tensor Hab reduces to the second order Lovelock tensor 

Rab = 9abR%B + (^RRab ~ ‘^RacbdR^^ + RacdeRb ~ “^Ra ^Rb^ 

= 2a {lab - ^Jab + Kab) , (A.4) 


where the quantities Rb, Jab, Kab are defined below. The linearized form of the curvatures are 

(i) Variation of Gab 

^Gab ^Rab 2 9ab ^R 2 ^ibR • (A.5) 

(ii) Variation of Rb (= R {Rab - i gabR)) 

SRb = Gab 5R + R 5Rab — ^ habR^ ■ (A.6) 

(hi) Variation of Jab (= RacbdR!"^ - i gabRcdR'^) 

SJab = SRapb^^RPq + RapbgSRP^ -^gabR^RRpq 

- [Rapb ''Rqr " ^gabRprRg " \ habRpqR^^ ■ (A.7) 

(iv) Variation of Kab (= RacdeRb - \ gabRcdefR^''^''^ - “^Ra "Rbc + 2RacbdR'"^) 

6Kab = 2R^a^1\,\SRt,)pg^ -^gabRP1R5Rpqr^ 

+ 2R^ dSRacb + ^RacbdSR^'^ - 4i?(, '^5R,)a 

— (^RacpqRbd^'^ — ^ dabRcpqrRd^'^^ + “^Racb ^Rde “ ‘^RacRbd^ 

- ^ habRpqrsR^^^^ • (A.8) 


The quantities defined with bar are to be taken about their background values. In the following, we 
adopt a slightly different scheme of linearization, which we find more convenient to use for a warped 
bulk geometry. The linearized equations take the form 


^Gab T bvf) SHab b , 


(A.9) 


where 


^Gab 


5Gab - GacK b = SGab - ^toSTab 


5Rab - - gabSR - RacK b 


(A.IO) 
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Slab = 5Iab-iach% = Gab5R + R5Rab-RRach‘^b (A.ll) 

SJab = SJab - Jach^ b = SRapb ^R^ , + Raphq^R^^ - \gabR^'^SRpq 

- (Rapb ^Rgr " ^9abRprRg " RapcgR^’^h ^ (A.12) 

6Kab = 6Kab-Kach^b 

= 2ii(„ ^Rb)pg " - IgabR^'^^ sSRpgr " + 2R'^ d^Racb ^ 

+ 2Racbd5R^^ -^Ria"5R^)c 

- (^RacpqRbd^'^ — 2 gabRcpqrRd^'^^ + 2Racb '^Rde “ 2RacRbd^ h‘^‘^ 

- [RapqrRc + ^Rapcq^ - 2Ra ^Rep) h " (A.13) 

In using (A.9), the RS fine-tuned relations do not appear explicitly in the equations of motion. 

Appendix B: Useful Identities with Warped Space-Time Metrics 

Let us consider the function defined by A{z) = log(A;| 2 :| -|- 1), where k in the inverse AdS curvature 
scale, k = 1/L. Since A{z) is a function of \z\, we have <92^ = yl'fl 2 | 2 :|, where dz\z\ = 20(z) — 1, A' 

is derivative of A with respect to its argument \z\, ©(z) is the Heaviside function. In D space-time 

dimensions, the warp factor A(z), as a general solution of the Einstein field equations modified by a 
Gauss-Bonnet term, reads 

y4(2:) = log /Ci I Zi I. (B.l) 

For simplicity we assume that ki = k 2 = ■ ■ ■ = k. A straight forward simplifications would rise to give 
the following results 

A'2 = e-2A(z) ^2 ^ [d^.\z\f = k^ {D - 4) 

i 

A" = -e-2^(^)fe2(D-4) + e-^(^)2A:^(5(zi). 

i 

Some useful identities that hold among the cross terms are 


d^.d^^A • d^iAd^^A = 2k^ ^ 5{zi) - k^ {D - 4)^ . (B.4) 

i 

d,.d,,Ad^id^^A = -e-‘^^^^Uk^^5{zi)5{zj) + {A'f . (B.5) 

i¥=j 

dzjdz^A + d^.Adz^A = e"^ k ^ {dz.d^^\zi\) = 2k5{zi) S^jzk ■ (B.6) 

i 

dz^^Ad'^'^hpu = kY^sgn{zi)dzihp^ ■ (B.7) 

i 

Zi\)d^^d^^hp^ = 2Y^5{zi)dl.hp^. (B.8) 

i i 


(B.2) 

(B.3) 
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Appendix C: Linear Expansions with the Gauss-Bonnet Term 

In the gauge = 0 = the first order linearized equations for h^i, take the general form 


= SGuu - =( -I dl-I dl + 




fiX 


]h 


Zi I ' l-flU 1 


= 2{D — A)ae 


2A 


(D-5) 


d,^Ad^^A-d,^d^^A] {dl + dl) 


+d^^A(^{D - 3)d,,d^'A - ^^ d,,Ad^*A^d: 
+ {dzAA + 9z,^9,,A)d^^d^^ 


h 


(C.l) 


where as dehned above i,j = 1,2, ■ ■ ■, (D — 4) count the number of extra (transverse) coordinates. In 
using the metric solution A{z) = log ^\^i\ linearized fluctuations simplify to 


1 

kd 


D-4 


— ^4 — + (L> — 2)ke ^ sgn{zi)d, 

+2a{D-A) 


h 


fiiy 


D-4 D-4 

{D - 3){D - 4)k^ (04 + □.) - 5{zi)n^ - 4ke^G) ^ sizi)dl 

i=l i^j 

D-4 D-4 1 

+{D - 3)^ (4 ^ 5{zi) -{D-2){D- 4)A:e-^(^)) sgn{z,)d,. V = > (C.2) 

i=l i=l ^ 

where = 9^^ +<9^2 ''' + ^zd- 4 ‘ redefining metric fluctuations as h^j,i, = one can 

remove from the hrst square bracket in Eq. (C.2) the single (linear) derivative term, and the kinetic 
term will have a canonical form. After this rescaling, linearized equations take the form 

\ r _ ! r> A\(r> ‘^\riu‘i 

Kd 

+(I)-4) 


-□4 -□„-(£)- 


{D - 2)ke-^G) ^ + {D-A){D -2)Dk 


h 


fAiy 


2{D -3){D- A)k‘^e-^^G) (n^ + n^) _ ^ 6{zi) ^4 

^ i 

^ 6{zi)d‘^^ + 8 (^{D - 3) ^ sgn{zi)5{zi)dzi - ^ 6{zj)sgn{zi)dzi 

i^j i i^j 

-4(E) - 2) ^ E>(19 - 3)(I) - 4)^ k^e-'^^G) + 2 ^ 5{zi) 5{zj) 


T _ 2 rp{m) 


- (E)2 -3D + 18) ^ 5{z 

i 

We have implemented these expressions in the bulk part of the text 


(C.3) 


Appendix D: Green Function in D Dimensions 

In Sec. 3, we made use of several analytic properties of bessel functions. Here we summarize some 
of the technical derivations that were involved. By defining in Eq. (23) Gp = (zz7E2)(^-b/2 and 
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, we arrive to 

(1 - 7 ) + zdz + - ■^— Qp{z,z') = Lz5{z - z'). (D.l) 

One can also derive the D dimensional boundary condition at z = L, analogous to the Neumann 
condition on the gravitational field that dzGD{x,x')\z=L = 0, to be 

q^z'^'^Gpiz, z%=L = 0 ■ (D.2) 

This is a boundary condition implied by continuity of the derivative fields on the brane. Equation (D.l) 
further implies the following matching conditions at z = z': 

G<\z=z' = G>\z=z' , dz{G> - Gc)\z=z' = - 7 ) (D.3) 

One can find the general solution of (D.l) by satisfying (D.3) and (D.2). The most general solution 
of (D.l), for z < z' \s 

gz<z' = iA{z') [(J ,_1 + X iqL)J,)Hi^\qz) - {Hj^l,{qL) + y (qL) Hi^^)Mqz)] , (D.4) 

where y = 7 /((z^ — 1) (1 — 7 )), = (D — l)/2, and Hl^\qL) = J^{qL) + iYy{qL) is the Hankel function 
of the first kind. For z > z', one can use the following boundary condition at the AdS horizon, z = 00 : 

Gz>z' = B{z!)BA^{qz). (D.5) 


Following [7, 9], we solve Eqs. (D.3,D.4,D.5)) and arrive to a D-dimensional Neumann propagator 
Qd(x,z\ x,J) = J^{qz^)H‘'J'>(qzy ) 


I \jz{az^)H7(qzy) 

( X-yqL) + x[qL)J iqL) \ 7 , 

\K\{qL)+xqLH7(qL)l J 


This gives the results in Ref. [9] for D = 5, and that of Ref. [7] with 7 = 0 and D = d -|- 1. When an 
argument of the propagator is at z' = L, (D. 6 ) can be reduced to Eq. (24) given in the text. When 
both arguments are on the brane, Eq.(24), after Bessel expansions and some algebraic manipulations, 
would simplify to yield 


Gnix. L; x',L)zz{l- J 1 


1-7 (P-3) [(P-4) + (P-6)7] 

1 -|- 7 qL 2{D — 5) 


X 


h^<,L + 0((,P)=)+...+ 


(1 - (qLf-' 

(1+7)2 Cl 



(D.7) 


where Ci is a dimensional constant (Ci = 1,4, 64, •••, for D = 5,7, 9, •••). The the last expression 
involving logarithmic term would be absent when the number of total dimensions D is even. 
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